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Using renormalization group methods we study multifractality in directed percolation. Our ap- 
proach is based on random lattice networks consisting of resistor like and diode like bonds with 
microscopic noise. These random resistor diode networks capture the features of isotropic as well 
as directed percolation. In this note we introduce a field theoretic Hamiltonian for the multifrac- 
tal properties at the transition from the non-percolating to the directed percolating phase. We 
investigate the multifractal moments of the current distribution and determine a family of critical 
exponents for these moments to two-loop order. 
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Nature is replete with systems where the exhaustive 
characterization of the distribution of a local physical 
quantity requires the introduction of an infinite set of 
independent critical exponents. For these systems the 
concept of multifractality H has been developed. Ex- 
amples for multifractal j^ystems include turbulence 
diffusion near fractals H, electrons in disordered me- 
dia polymers in disordered media random fer- 
romagnets ||, chaotic dissipative systems 0, random 
resistor networks §, § [k], [ll} 0, [jj) , and heartbeat @. 

Here, we discuss multifractality in the context of di- 
rected percolation (DP) |^5[ . Consider two terminal 
points x and x' on a fractal DP cluster. A current / 
is inserted at x and withdrawn at x' . Particularly inter- 
esting quantities are the moments 
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of the currents J& in all the bonds b on the cluster. Specif- 
ically, Mj°\ Mj, Mj 2 \ and m] 00 -* are proportional to 
the number of bonds on the backbone (bonds which carry 
nonzero current), the total resistance, the noise (second 
cumulant of the resistance fluctuations), and the number 
of the red bonds (bonds which carry the full current). 
If one measures Mj ' in the direction along which the 
cluster is oriented, one finds that 
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where — xi | is the distance between the two terminals 
and v\\ is the longitudinal correlation length exponent of 
the DP universality class. The multifractality manifests 
itself in the nonlinear dependence of tpi on I. This behav- 
ior is much richer and more interesting than the usual gap 
scaling commonly encountered in critical phenomena. Of 
course, it is a challenging task to characterize the current 
distribution because one has to determine a multitude of 
critical exponents ipi. 



In this paper we study multifractality in DP based on 
the random resistor diode network (RDN) introduce d by 
Redner H). Guided by ideas of Stephen JM, Harris @, 
and Park, Harris and Lubensky (PHL) ]10| we introduce 
a field theoretic Hamiltonian for RDN with noise. We de- 
termine the multifractal exponents ipi by computing the 
scaling indices of noise cummulants Cr which are known 

to be proportional (cf. Ref. to the Mf 1 by virtue of 
Cohn's theorem pQ]. A central role in this analysis is 
played by our concept of master operators jO, 12, b|. 

A RDN consists of a d-dimensional hypercubic lattice 
in which nearest-neighbor sites are connected by a resis- 
tor, a positive diode (conducting only in a distinguished 
direction), a negative diode (conducting only opposite to 
the distinguished direction), or an insulator with respec- 
tive probabilities p, p + , p_ , and q = 1 — p — p+ — p~. 
In the three dimensional phase diagram (pictured as a 
tetrahedron spanned by the four probabilities) one finds 
a nonpercolating and three percolating phases The 
percolating phases are isotropic, positively directed, or 
negatively directed. Between the phases there are sur- 
faces of continuous transitions. All four phases meet 
along a multicritical line, where < r := p + = p- < 1/2 
and p = p c (r). On the entire multicritical line, i.e., inde- 
pendently of r, one finds the scaling properties of usual 
isotropic percolation (r = 0). For the crossover from 
IP to DP see, e.g., Ref. |21|. In this paper we focus on 



the vicinity of the critical surface separating the non- 
percolating and the positively directed phase. Here, typ- 
ical clusters are anisotropic and they are characterized 
by two different correlation lengths: £|| (parallel to the 
distinguished direction) and £± (perpendicular to it). As 
one approaches the critical surface, the two correlation 
lengths diverge with the exponents and v± of the DP 
universality class. 

We choose n = l/yd(l, . . . , 1) for the distinguished 
direction, and we assume that the bonds j\ between 
two nearest neighboring sites i and j are directed so 
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n > 0. The directed bonds obey the non- 
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linear Ohm's law erj, (V&) Vb = h, where Vb is the po- 
tential drop over the bond and lb denotes the current 
flowing from j to i. The bond conductances Gb = ft 76 
are equally and independently distributed random vari- 
ables. The 76 take on the values 1, 6 (V), 9 (—V), and 
with respective probabilities p, p+, p_, and q. 9 denotes 
the Heaviside function. Static noise is introduced by 
distributing the q, according to some distribution func- 
tion / with mean ? and higher cumulants A''- 2 ) sat- 
isfying AW <C S l ■ The noise average will be denoted 
by {•••}/ = / lib fob f (o>) ■ ■ ■ and its lih cumulant by 

{• • •'}y C ' 1 . Both kinds of disorder, the random dilution of 
the lattice and the noise, influence the statistical proper- 
ties of the resistance R+(x, x') [p2| between two sites x 
and x' . They are reflected by the noise cumulants 

Cr{x,x) = — J - , 3 

where ( • • • ) c denotes the average over the configurations 
C of the randomly occupied bonds and x+{x,x') is an 
indicator function which is unity if x and x' are positively 
connected (i.e., if current can percolate from a; to a;') and 
zero otherwise. 

Now we determine ipi by exploiting the proportionality 
of Mj and . Our actual strategy is the following: (i) 

we devise a generating function for the , (ii) we calcu- 
late this generating function by renormalized field theory, 

and (iii) we extract (7^ from the generating function. 

To facilitate the averaging procedure we employ the 
replica technique Following PHL we replicate 

the voltages DxE-fold, V x -> V x = (Vi (Q ' ,3) )f; / f =1 . 
In order to set up the desired generating function we 

introduce ip^(x) — exp(i\ ■ V x ), where A • V x — 
Ea'^=i XS a '^V x a '^ and A ^ and define corresponding 
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G(x,x',X) = 
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Here P({V}) = E^ =1 E^^H^ is the 
replicated power dissipated on the backbone and Z is the 
usual normalization. 

Due to the appearance of the ^-functions in the power 
P the integration in Eq. (|J) is not Gaussian. However, 
it can be carried out in an approximating manner by 
applying the saddle point method as it was done by 
Harris jl8j in the related problem of nonlinear random 
resistor networks. We extract the leading contribution 
to the integral steaming from the maximum of the in- 
tegrand. This maximum is determined by the solution 
of Kirchhoff 's equations, i.e., by the resistance R+(x,x'). 



Hence one finds (details on the steps can be gleaned from 
Refs. H H)) that 
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where K t (\) = Ef=i[E„=i(^ /3) ) 2 ] i '■ Equation (§) 
shows that G represents the desired generating function 
from which C^can be calculated via taking the deriva- 
tive with respect to Ki(\). 

A glance at Eq. (0) reveals the benefit of the replica- 
tion. It provides an additional parameter D which we 
can tune to zero. In this replica limit the normalization 
denominator in Eq. (Q) is taken to unity and in particular 
it no longer depends on the disorder. Hence we are left 
with the much simpler task to average the exponential in 
Eq. (0) which gives us an effective Hamiltonian 
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For technical reasons 10, [llj |l2|,Jl3[ we switch to dis- 

cretized voltages "& and currents A taking values on a 
discrete D x -©-dimensional torus. For the saddle point 
method to be reliable we work near the limit when all the 
components of A are equal and continue to large imagi- 
nary values. Accordingly we set E3| A^'' 3 ' = i\q + £( q -' 3 ) 
with real Ao and ^ a -^\ Ea 'JLi ^"'^ = 0, and impose 

the conditions A 2 <C D~ x and £ 2 <C 1. 

To refine H rcp towards a field theoretic Hamiltonian, 
we expand -ff rop in terms of (x). The steps are analo- 
gous to those in Ref. |TJ| and are skipped here for brief- 
ness. The so obtained expression is converted into a 
Landau-Ginzburg- Wilson-type functional 
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by applying the usual coarse graining procedure. The pa- 
rameter t specifies the "distance" from the critical sur- 
face under consideration. The vector v lies in the distin- 
guished direction, v = vn. r and v depend on the three 
probabilities p, p+, and p_. w and vi are the coarse 
grained analogs of q^ 1 and A^/s 21 , respectively. In the 
limit w — > and m — > our Hamiltonian TL describes the 
usual purely geometric DP. Indeed Ti leads for w — > 
and v\ — > to exactly the same perturbation series as 
obtained in Refs. 121, M, E^l. 
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Now we address the relevance of the vi . A straightfor- 
ward scaling analysis reveals that 



cf((x,x');T,v, 
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where /; is a scaling function. Note that the coupling 
constants Vk appear only as Vk/w k . Dimensional analysis 

shows that wX 2 ~ p? and VkKk(X) ~ /J- 2 , where [i is the 
usual inverse length scale. Thus, Vk/w k ~ /i 2 ~ 2fe , i.e., the 
Vk/w have a negative naive dimension which decreases 
drastically with increasing k. We conclude that the Vk 
are irrelevant couplings. Though irrelevant, we must not 
set Vk = in calculating the multifractal exponents. To 
see this we expand the scaling function /; yielding 



C«((x,x');t,v, 
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with C^' being an expansion coefficient depending on x, 

x', r, and v. Note that = because the correspond- 
ing terms are not generated in the perturbation calcula- 
tion. The first term on the right hand side of Eq. (||) gives 

the leading behavior. Thus, vanishes upon setting 
vi = and we cannot gain any further information about 
C}j . In other words, the vi are dangerously irrelevant in 

investigating the critical properties of the C^~ 2 ' • 

We apply standard methods of field theory Q and 
perform a diagrammatic perturbation calculation up to 
two-loop order. Since the coupling constants vi are ir- 
relevant, they cannot be treated in the same fashion as 
the other coupling constants r and w pertaining to the 

bilinear part of H. Using [r + p 2 + wX 2 + J^i viKi(X)]^ 1 
as the Gaussian propagator would poison the Feynman 
diagrams. This can be understood by expanding any 
given diagram in terms of vi. For increasing orders of 
this expansion one encounters increasing orders of primi- 
tive divergence. Hence one has to truncate this expansion 
at linear order in vi which is equivalent to treating vi by 
means of the insertion 



0(9 



1 



^Ei^(^i(p)^(-p), (io) 



where </>^(p) denotes the Fourier transform of ip<?(x). 
This insertion does not only generate primitive divergen- 
cies proportional to A/ (A), but also primitive divergences 
corresponding to all operators Of of the generic form 

X 2a p 2b (f> n having the same or a lower naive dimension 
than 0(9. In other words: all operators of the form 

X 2a p 2b (j) n with a + b + n < 1 + 2 are generated. Those 
0( associated with a lower naive dimension than 0(9 
can be neglected because they just lead to corrections 
which vanish at the critical point (cf. Ref. p6| ). Anyway, 
one still needs a myriad of operators as counterterms in 



the Hamiltonian. Inserting either of these 0(9 , however, 
does not generate 0(9 because the O^p possess of higher 
symmetries than 0(9. Thus, the renormalization scheme 
is given by 

OIL =Z®0® +£^ )0 ?rc„, (11) 
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and one solely needs in calculating the scaling index 
of r ( 9 n . What we encounter here is another realization 
of our concept of master operators (llj [l^, O. Each 

multifractal moment Mj corresponds one-to-one to a 
dangerously irrelevant master operator 0(9. Though a 
myriad of servant operators is involved in the renormal- 
ization of the masters 0(9 the scaling behavior of Mj is 
governed by 0(9 only. The servant operators can be ne- 
glected in determining the scaling index of their master 
operator, i.e., we are spared the computation and diago- 
nalization of giant renormalization matrices. 

By employing dimensional regularization and mini- 
mal subtraction we proceed with standard techniques of 
renormalized field theory. We calculate to two-loop 
order. From the renormalization group equation and di- 
mensional analysis we deduce that the correlation func- 
tion G scales at criticality as 



G (0,x, A J = x 



7 



1 + WX Xn 



x/.^+^p)^/,^)}^) 

where higher order terms have been discarded, rj, m , 
and z = Vn/fx are the critical exponents for DP known 
to second order in e = 5 — d ]24|, [25|. f,f w , and /„ are 
scaling functions. 4> is t ne resistance exponent given to 
second order in e in Ref. |27[ . In Eq. ([l3l) we introduced 
the multifractal exponents ipi = 




where 
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stands for the Wilson function 7(9 — 



-M^lnZ(9| 



evaluated at the infrared stable fixed point. By taking 
the derivative with respect to A"; (A) we can now de- 
duce the scaling behavior of the multifractal moments: 

if measured parallel to the distinguished 
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direction. For measurements in other directions it is ap- 
propriate to choose a length scale L and to express the 
longitudinal and the transverse coordinates in terms of 
L: Ixi I ~ L and xn ~ L z . With this choice we find that 
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The result for the multifractal exponents remains to 
be stated. For I > we obtain to second order in e 



xpi = l + 
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with a(l) and b(l) taking on the values listed in Table 1. 
tpo and ipi stem from extending the sum over I in Eq. (pf) 



4 



to comprise I = and 1 = 1. Equation ( |l4| ) fulfills several 
consistency checks, i/jq is related to the fractal dimension 
Db of DP clusters via Db = l+ipo/u±. Equation (|l4|) is 
in agreement with the e-expansions of u± [§3 ^| and D b 
(cf. Ref. |2?| ) to second order in e. ipi is in conformity 
with our result for the resistance exponent <f> given in 
Ref. ^8) . As anticipated ipi is for reasonable values of 
e a convex monotonically decreasing function of I |p9| . It 
tends to unity for large I as one expects from the relation 
of ipoo to the fractal dimension of the singly connected 
(red) bonds, cf. Refs. §, [3§. 

While directed percolation is a vivid area of research 
since about 20 years and multifractality in ordinary per- 
colation is studied since about 15 years, the problem of 
multifractality in DP has not been addressed hitherto, 
at least to our knowledge. In the present note, we intro- 
duced a field theoretic Hamiltonian for the multifractal 



properties of resistor diode percolation. We applied our 
concept of master operators. It turned out that it works 
consistently as a tool to describe the multifractal prop- 
erties of DP by renormalized field theory. Without the 
concept of master operators which avoids the calculation 
and diagonalization of giant renormalization matrices the 
actual calculations are hardly feasible. We introduced 
the family of multifractal exponents {ipi} in the spirit of 
PHL. We calculated ipi for non-negative I to two-loop or- 
der. Our result fulfills several consistency checks. It is 
certainly desirable to have numerical data for compari- 
son with our analytic results, in particular in dimensions 
to which the e-expansion can be continued reliably. We 
hope that this paper triggers such simulations. 
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Deutsche Forschungsgemeinschaft. 
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TABLE I: The coefficients a(l) and b(l) appearing in Eq. 
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